Method of transmitting data using constellation rearrangement

ABSTRACT

A data retransmission method using hybrid automatic repeat request (HARQ) includes transmitting a data block, receiving a retransmission request signal for the data block, generating a retransmission block by performing swapping or inversion between bits constituting the data block according to the retransmission request signal, and transmitting the retransmission block.

TECHNICAL FIELD

The present invention relates to wireless communications, and more particularly, to a method and apparatus using constellation rearrangement in a wireless communication system.

BACKGROUND ART

Current development in advanced wireless communication has led to the requirement of high spectral efficiency and reliable communication. Unfortunately, packet errors by fading channel environment and interferences originated from various sources make the capacity of overall system to be limited.

Hybrid Automatic Repeat Request (HARQ) which is ARQ protocol combined with Forward Error Correction (FEC) is strongly considered as one of cutting edge technologies for future reliable communication. The HARQ scheme can largely be classified into the type of two. One is HARQ-Chase Combining (CC) which is disclosed in D. Chase, Code Combining: A maximum-likelihood decoding approach for combining an arbitrary number of noisy packets, IEEE Trans. on Commun., Vol. 33, pp. 593-607, May 1985. The other scheme is HARQ-Increment Redundancy (IR). In the HARQ-CC, when a receiver detects an error through cyclic redundancy checking (CRC) while decoding the transmitted packet, the same packet with the same modulation and coding is sent to the receiver, repeatedly. Meanwhile, HARQ-IR retransmits different packets in order to achieve the coding gain, in which parity bits can be manipulated through puncturing and repetition.

Multiple Input Multiple Output (MIMO) systems are regarded as one of the most promising research areas of wireless communication. Spatial diversity provided by multiple antenna configurations for both transmitter and receiver is known to tremendously increase system capacity without additional bandwidth. As a result, various approaches have been studied to use the benefit of transmit diversity and received diversity.

Constellation Rearrangement disclosed by the PCT international application No. PCT/KR2007/003625 filed by this applicant provides additional gain through averaging the difference of inherent reliability between component bits.

There is a need to improve performance of system by applying constellation rearrangement to various schemes such as MIMO system, HARQ, multi-level modulation, etc.

DISCLOSURE OF INVENTION Technical Problem

The present invention provides a method and apparatus for performing constellation rearrangement in a wireless communication system.

The present invention also provides a method of performing constellation rearrangement in a multiple-antenna system using a multi-modulation scheme.

Technical Solution

In an aspect, a data retransmission method using hybrid automatic repeat request (HARQ) includes transmitting a data block, receiving a retransmission request signal for the data block, generating a retransmission block by performing swapping or inversion between bits constituting the data block according to the retransmission request signal, and transmitting the retransmission block.

The swapping between bits may be performed by swapping a bit having a high bit reliability and a bit having a low bit reliability in previous transmission.

In addition, the data block may be a transmission block having a matrix format and consisting of rows whose number is equal to the number of transmit antennas, each row of the transmission block may be transmitted through one transmit antenna, bits constituting one row may represent bits of one data symbol, and the data symbol may be modulated using at least two modulation schemes.

In another aspect, a data retransmission method in a multiple-antenna system includes transmitting a first data symbol modulated with a first modulation scheme through a first transmit antenna and transmitting a second data symbol modulated with a second modulation scheme through a second transmit antenna, generating a new first data symbol and a new second data symbol by performing swapping or inversion on bits constituting the first data symbol and bits constituting the second data symbol, and transmitting the new first data symbol through the first transmit antenna and transmitting the new second data symbol through the second transmit antenna.

In still another aspect, an apparatus for wireless communication includes an radio frequency (RF) unit for transmitting and receiving a radio signal, and a processor coupled with the RF unit, for performing HARQ-increment redundancy (IR), and configured to sequentially configure data blocks in a cyclic buffer of the HARQ-IR according to a retransmission request, and transmit the data blocks, wherein if all or some parts of the data block are wrapped around with a previously transmitted data block, constellation rearrangement is performed on the wraparound parts.

ADVANTAGEOUS EFFECTS

An additional diversity gain can be obtained in a multiple input multiple output (MIMO) system by using constellation rearrangement. Therefore, performance of a wireless communication system can be improved.

BRIEF DESCRIPTION OF THE DRAWINGS

FIG. 1 is a block diagram showing a transmitter and a receiver according to an embodiment of the present invention.

FIG. 2 shows an example of an operation of an adaptive mapper for multi-modulation transmission in case of using two transmit antennas.

FIG. 3 is a flow diagram showing a data transmission method using a wireless communication system of FIG. 1.

FIG. 4 shows a transmission block in case of using multiple antennas.

FIG. 5 shows a signal constellation of an M-ary quadrature amplitude modulation (M-QAM) scheme.

FIG. 6 is a graph showing changes in a minimum squared error (MSE) of average bit reliability with respect to the number of retransmissions according to a proposed bit swapping scheme.

FIG. 7 shows a bit swapping and inversion (BSI) scheme in a 2×2 multiple input multiple output (MIMO) system using 16-QAM and 64-QAM according to an embodiment of the present invention.

FIG. 8 shows an example of a bit shuffling between antennas (BSA) scheme in single modulation transmission.

FIG. 9 shows a BSA scheme in a 5×5 MIMO system using 16-QAM and 64-QAM according to an embodiment of the present invention.

FIG. 10 shows an example of applying a BSI scheme and a BSA scheme to a 2×2 MIMO system using 16-QAM and 64-QAM.

FIG. 11 is a graph showing a result obtained by performing simulations in a 2×2 MIMO system using 16-QAM and 64-QAM.

FIG. 12 shows gray mapping for 16-QAM.

FIG. 13 shows decision boundary to calculate a bit error probability.

FIG. 14 shows a swapping operation.

FIG. 15 shows an inversion operation.

FIG. 16 shows decision boundary of 8-pulse amplitude modulation (PAM).

FIG. 17 shows frame error rate (FER) performances according to the number of bits shuffled between transmit antennas.

FIG. 18 shows a comparison result of FER performances of optimal bitwise mapping and conventional mapping.

FIG. 19 shows hybrid automatic repeat request (HARQ)-incremental redundancy (IR) using a swapping operation.

FIG. 20 shows HARQ-IR using a swapping operation and an inversion operation.

FIG. 21 is a block diagram showing an apparatus for wireless communication according to an embodiment of the present invention.

MODE FOR THE INVENTION

The technology described below can be used in various wireless communication systems such as code division multiple access (CDMA), frequency division multiple access (FDMA), time division multiple access (TDMA), orthogonal frequency division multiple access (OFDMA), single carrier frequency division multiple access (SC-FDMA), etc. The CDMA can be implemented with a radio technology such as universal terrestrial radio access (UTRA) or CDMA-2000. The TDMA can be implemented with a radio technology such as global system for mobile communications (GSM)/general packet ratio service (GPRS)/enhanced data rate for GSM evolution (EDGE). The OFDMA can be implemented with a radio technology such as institute of electrical and electronics engineers (IEEE) 802.11 (Wi-Fi), IEEE 802.16 (WiMAX), IEEE 802-20, evolved UTRA (E-UTRA), etc. The UTRA is a part of a universal mobile telecommunication system (UMTS). 3rd generation partnership project (3GPP) long term evolution (LTE) is a part of an evolved UMTS (E-UMTS) using the E-UTRA. The 3GPP LTE uses the OFDMA in downlink and uses the SC-FDMA in uplink. LTE-advance (LTE-A) is an evolution of the 3GPP LTE.

This technology can be used in downlink or uplink. In general, a downlink denotes a communication link from a base station (BS) to a user equipment (UE), and an uplink denotes a communication link from the UE to the BS. The BS is generally a fixed station that communicates with the UE and may be referred to as another terminology, such as a node-B, a base transceiver system (BTS), an access point, etc. The UE may be fixed or mobile, and may be referred to as another terminology, such as a mobile station (MS), a user terminal (UT), a subscriber station (SS), a wireless device, etc.

The communication system may be a multiple-antenna system having a plurality of transmit antennas. Hereinafter, a multiple-input multiple-output (MIMO) system denotes a system using a plurality of transmit antennas and/or a plurality of receive antennas.

FIG. 1 is a block diagram showing a transmitter and a receiver according to an embodiment of the present invention. A transmitter 100 and a receiver 200 implement hybrid automatic repeat request (HARQ). The transmitter 100 and the receiver 200 can be regarded as a transceiver that performs both a transmission function and a reception function. For clear explanation of data retransmission, one side that transmits and retransmits data is referred to as the transmitter 100, and the other side that receives data and requests data retransmission is referred to as the receiver 200. In downlink, the transmitter 100 may be a part of a BS, and the receiver 200 may be a part of a UE. In uplink, the transmitter 100 may be a part of the UE, and the receiver 200 may be a part of the BS. The BS may include a plurality of receivers and a plurality of transmitters. The UE may include a plurality of receivers and a plurality of transmitters.

Referring to FIG. 1, a transmitter 100 includes a channel encoder 110, an adaptive mapper 120, a spatial encoder 130, a controller 150, and a receive circuitory 180. Further, the transmitter 100 includes Nt modulators 140-1, . . . , 140-Nt and Nt transmit antennas 190-1, . . . , 190-Nt, where Nt is greater than one (i.e., Nt>1).

The channel encoder 110 receives a stream of information bits and encodes the received stream of information bits according to a predetermined coding scheme. As a result, coded data is generated. The adaptive mapper 120 modulates the coded data according to a predetermined modulation scheme and thus provides a data symbol. The adaptive mapper 120 can use at least two modulation schemes. The adaptive mapper 120 maps the coded data to the data symbol representing a position on a signal constellation. Further, the adaptive mapper 120 adaptively remaps the coded data in response to a retransmission request message of the controller 150. There is no limit in the modulation scheme used by the adaptive mapper 120. The modulation scheme may be an M-ary quadrature amplitude modulation (M-QAM). Examples of the M-QAM include 16-QAM, 64-QAM, and 256-QAM. Detailed operations of the adaptive mapper 120 will be described below.

The spatial encoder 130 processes data symbols output through the adaptive mapper 120 according to a MIMO pre-processing scheme. The modulators 140-1, . . . , 140-Nt modulate symbols output from the spatial encoder 130 and transmit the modulated symbols through the respective transmit antennas 190-1, . . . , 190-Nt. When the modulators 140-1, . . . , 140-Nt perform an inverse fast Fourier transform (IFFT), orthogonal frequency division multiplexing (OFDM) symbols are output. The receive circuitory 180 receives signals transmitted from the receiver 200 through the transmit antennas 190-1, . . . , 190-Nt. The receive circuitory 180 digitizes the received signals and then transmits the digitized signals to the controller 150.

The controller 150 controls overall operations of the transmitter 100. The controller 150 extracts information from signals received from the receive circuitory 180. An operation of extracting the information includes general modulation and decoding. The extracted information may include a retransmission request signal. The controller 150 prepares a retransmission symbol by controlling the adaptive mapper 120 in response to the retransmission request signal.

A channel quality indicator (CQI) may be included in the information extracted from the signal received from the receive circuitory 180. The CQI may be information on a channel condition from the receiver 200 to the transmitter 100 or index information on a modulation and coding scheme. The CQI can be used by the controller 150 to control the channel encoder 110 or the adaptive mapper 120. Thus, a coding scheme of the channel encoder 110 or a mapping scheme of the adaptive mapper 120 can be adaptively changed.

Meanwhile, a receiver 200 includes a spatial decoder 220, a demapper 230, a channel decoder 250, an error detector 260, a controller 270, and a transmit circuitory 280. Further, the receiver 200 includes Nr receive antennas 290-1, . . . , 290-Nr, where Nr is greater than one (i.e., Nr>1).

Signals received from the receive antennas 290-1, . . . , 290-Nr are demodulated by demodulators 210-1, . . . , 210-Nr and then are input to the spatial decoder 220. The spatial decoder 220 processes the received signals according to a MIMO post-processing scheme in response to a MIMO control signal. The MIMO control signal controls decoding according to a space time coding (STC) scheme of the receiver 200. The MIMO control signal may be pre-defined by a memory (not shown) of the controller 270. Alternatively, the MIMO control signal may be received from the transmitter 100.

The demapper 230 demaps data symbols from coded data according to a demapping control signal provided from the controller 270. The demapping control signal controls the demapper 230 according to a mapping scheme used in the adaptive mapper 120 of the transmitter 100. The demapping control signal may be pre-defined by the memory of the controller 270. Alternatively, the demapping control signal may be received from the transmitter 100.

The receiver 200 may include a combination unit 240 that combines a retransmitted symbol and a previous symbol. That is, in case of using an HARQ scheme such as HARQ-chase combining (CC) or HARQ-incremental redundancy (IR), the combination unit 240 combines retransmitted symbols and previous symbols. A combining scheme in use may be an equal-gain combining scheme in which combination is performed using an average value by assigning the same weight factor to both previous data and retransmitted data. The combining scheme may be a maximal ratio combining (MRC) scheme in which weight factors are assigned to respective pieces of data. There is no limit in the combining scheme, and thus other various schemes may also be used.

The technical features of the present invention are not limited to the HARQ-CC or HARQ-UR scheme, and can also apply to an HARQ scheme in which channel decoding is performed using only retransmitted symbols without being combined with previous symbols. In this case, the receiver 200 may not include the combination unit 240 as indicated by a dotted line in the figure.

The channel decoder 250 decodes coded data according to a predetermined decoding scheme. The error detector 260 detects an error from a decoded data bit by using cyclic redundancy check (CRC).

The controller 270 controls overall operations of the receiver 200 and provides a retransmission request signal or the like to the transmit circuitory 280. For this, the controller 270 can perform general channel encoding, modulation, etc. The controller 270 receives a result of error detection from the channel decoder 250 and determines whether to request retransmission. The controller 270 may generate a positive acknowledgement (ACK) signal if no error is detected, and may generate a negative acknowledgement (NACK) signal if an error is detected. The NACK signal may be the retransmission request signal.

Further, the controller 270 can provide a CQI signal by measuring channel quality from received signals. The CQI signal is a feedback signal to be fed back to the transmitter 100. The feedback signal indicates channel quality such as a signal-to-noise ratio (SNR) or an error rate. The transmit circuitory 280 receives the retransmission request signal or the like from the controller 270, and transmits the received signal through the receive antennas 290-1, . . . , 290-Nr.

I. Constellation Rearrangement for Multi-Modulation Transmission.

Different modulation schemes (hereinafter, multi-modulation transmission) can be used for respective transmit antennas in a MIMO system. For example, a 1st transmit antenna may use 16-QAM and a 2nd transmit antenna may use 64-QAM. A method of performing constellation rearrangement in a multiple-antenna system in which various modulation schemes co-exist for respective transmit antennas will be described.

FIG. 2 shows an example of an operation of an adaptive mapper for multi-modulation transmission in case of using two transmit antennas.

Referring to FIG. 2, the 1st transmit antenna 190-1 uses 16-QAM and the 2nd transmit antenna 190-2 uses 64-QAM. In initial transmission, four bits constitute one data symbol in the 1st transmit antenna 190-1, and six bits constitute one data symbol in the 2nd transmit antenna 190-2. Data symbols for initial transmission in the 1st transmit antenna 190-1 are expressed by {b₁, b₂, b₃, b₄}. Data symbols for initial transmission in the 2nd transmit antenna 190-2 are expressed by {b₅,b₆,b₇,b₈,b₉,b₁₀}. When an error occurs according to a channel condition, data symbols transmitted from the respective transmit antennas are retransmitted. The adaptive mapper 120 configures bits constituting a retransmission symbol for each antenna by performing swapping or inversion of each bit on a signal constellation of the symbol in comparison with initial transmission. Unlike the initial transmission, in retransmission, the 1st transmit antenna 190-1 uses 64-QAM and the 2nd transmit antenna 190-2 uses 16-QAM.

According to this example, in retransmission, data symbols of the 1st transmit antenna 190-1 are remapped to {b₃, b ₅, b₁, b ₇}, and data symbols of the 2nd transmit inversion operation.

FIG. 3 is a flow diagram showing a data transmission method using the wireless communication system of FIG. 1.

Referring to FIG. 3, the transmitter 100 transmits a transmission (Tx) block S^((O)) (step S110). The Tx block is formed in an Nt×2B_(max) matrix format, and is a data block max mapped onto a signal constellation by the adaptive mapper 120. Nt denotes the number of transmit antennas. 2B_(max) denotes the number of index bits for a modulation scheme having a highest order among modulation schemes used in multi-modulation transmission. If a maximum modulation order is M-QAM, 2B_(max) is log₂M. For example, 2B_(max) is 4 in 16-QAM. The Tx block expressed in a mathematical form consists of a predetermined number of rows, wherein the predetermined number is equal to the number of transmit antennas. Each row is transmitted through one transmit antenna. Bits constituting one row represent bits of one data symbol. In one Tx block, the data symbol is modulated using at least two modulation schemes.

The superscript of the Tx block S denotes the number of retransmissions. For example, S⁽⁰⁾ denotes a Tx block for initial transmission, and S⁽¹⁾ denotes a Tx block for 1st retransmission.

The receiver 200 detects an error from the received Tx block S⁽⁰⁾ (step S120). If no error is detected, the receiver 200 transmits an ACK signal to the transmitter 100, and waits for transmission of a next Tx block. However, it will be assumed herein that the receiver 200 detects the error and thus transmits a NACK signal as a retransmission request signal (step S130).

Upon receiving the NACK signal, the transmitter 100 transmits a retransmitted Tx (ReTx) block S⁽¹⁾ (step S140). Upon receiving the NACK signal, the controller 150 controls the adaptive mapper 120 to remap the Tx block S₍₀₎ in a bitwise manner and/or in a spatial manner, and thus configures the ReTx block S⁽¹⁾. Various schemes can be used as a remapping scheme used in retransmission, which will be described below.

The receiver 200 detects an error from the received ReTx block S⁽¹⁾ (step S150). In this case, the combination unit 240 can combine the previous Tx block S⁽⁰⁾ and the ReTx block S⁽¹⁾.

If no error is detected, the receiver 200 transmits an ACK signal to the transmitter 100, and waits for transmission of a next symbol. However, it will be assumed herein that the receiver 200 detects the error and thus transmits a NACK signal as a retransmission request signal (step S160).

Upon receiving the NACK signal, the transmitter 100 transmits a remapped ReTx block S⁽²⁾ (step S170). The adaptive mapper 120 remaps the Tx block S⁽⁰⁾ in a bitwise manner and/or in a spatial manner, and thus configures the ReTx block S⁽²⁾.

The receiver 200 detects an error from the received ReTx block S⁽²⁾ (step S180). The receiver 200 transmits an ACK signal or a NACK signal to the transmitter 100 according to an error detection result (step S190). When the ACK signal is transmitted, retransmission for a corresponding Tx block is finished. A retransmission request in response to the NACK signal can be repeated up to a predetermined number M of times, where M is an iteration number greater than 0 (i.e., M≧1). If errors are continuously detected even after M-th retransmission is performed, a retransmission process can be reset and then a next Tx block can start to be transmitted. Alternatively, transmission of a current Tx block may be resumed.

FIG. 4 shows a Tx block in case of using multiple antennas.

Referring to FIG. 4, T₀ denotes initial transmission, T₁ denotes 2nd transmission, that is, 1st retransmission, and T_(m) denotes (m+1)th transmission, that is, m-th retransmission. S⁽⁰⁾ denotes a Tx block for initial transmission. S⁽¹⁾ denotes a Tx block for 1st retransmission. S^((m)) denotes a Tx block for m-th retransmission.

If different modulation schemes are used by the respective transmit antennas, and if 2B_(n) denotes the number of index bits for a modulation scheme of an n-th transmit antenna, then an n-th row s⁽⁰⁾ constituting the Tx block S⁽⁰⁾ can be expressed with bits indicating an in-phase (I)-axis and a quadrature (Q)-axis as shown:

MathFigure 1

s_(n) ⁽⁰⁾=[i_(n,1) . . . i_(n,B) _(n) q_(n,1) . . . q_(n,B) _(n) ]  [Math.1]

where i and q respectively denote bits indicating the I-axis and the Q-axis on the signal constellation. The positions of the I-axis and the Q-axis are not absolute positions. That is, if one axis on the signal constellation is referred to as the I-axis, the other axis is referred to as the Q-axis.

Thus, the Tx block S⁽⁰⁾ can be expressed as follows:

$\begin{matrix} {{MathFigure}{\mspace{11mu} \;}2} & \; \\ \begin{matrix} {S^{(0)} = \begin{bmatrix} I^{(0)} & Q^{(0)} \end{bmatrix}} \\ {= \begin{bmatrix} i_{1,1} & \ldots & i_{1,B_{1}} & x & x & q_{1,1} & \ldots & q_{1,B_{1}} & x & x \\ i_{2,1} & i_{2,2} & \ldots & i_{1,B_{2}} & x & q_{2,1} & q_{2,2} & \ldots & q_{1,B_{2}} & x \\ \vdots & \vdots & \vdots & \vdots & x & \vdots & \vdots & \vdots & \vdots & x \\ i_{{Nt},1} & i_{{Nt},2} & i_{{Nt},3} & \ldots & i_{{Nt},B_{\max}} & q_{{Nt},1} & q_{{Nt},2} & q_{{Nt},{3q}} & \ldots & q_{{Nt},B_{\max}} \end{bmatrix}} \end{matrix} & \left\lbrack {{Math}.\mspace{14mu} 2} \right\rbrack \end{matrix}$

where Nt denotes the number of transmit antennas, and x denotes an empty element which indicates that no value exists in a corresponding location. Likewise, retransmission blocks S⁽¹⁾, S^((m)) can be expressed with an Nt×2B_(max) matrix.

In a case where the respective transmit antennas use different modulation schemes or the same modulation scheme, the retransmission blocks S⁽¹⁾, S^((m)) are searched for, wherein bit swapping and inversion are performed on these blocks to reduce a bit error probability in each retransmission.

A link performance gain can be obtained in a retransmission method based on the adaptive mapper for the following two reasons. First, a position of each bit has an unequal bit importance due to a QAM characteristic. When diversity can be obtained by varying mapping of the signal constellation, the diversity is called mapping diversity. The mapping diversity is obtained by performing swapping of bits constituting a data symbol or by performing bit inversion when retransmission is made. Second, spatial diversity can be obtained by shuffling transmit antennas when retransmission is made using multiple antennas.

Bit swapping and inversion (BSI) denotes horizontal rearrangement for obtaining the mapping diversity. Bit shuffling between antennas (BSA) denotes vertical rearrangement for obtaining the spatial diversity. If all transmit antennas use the same modulation scheme, bit mappings obtained using the BSI and the BSA are independent from each other. On the other hand, if each transmit antenna uses a different modulation scheme, the obtained bit mappings are not completely independent since swapping between bits located in different rows has to be considered using the BSI.

There is a need for a mapping scheme capable of optimizing the bit error probability and the spatial diversity.

FIG. 5 shows a signal constellation of an M-QAM scheme. It is assumed that the location of signal constellation conforms to general gray mapping.

Referring to FIG. 5, if M_(n)-QAM denotes a modulation scheme used in an n-th transmit antenna, a total of M_(n) signal constellations exist, and the number of index bits for each modulation scheme is 2B_(n)=log₂M_(n).

From the perspective of the I-axis, if D_(n) denotes a minimum distance between positions on the signal constellation, a position I(c) of an arbitrary Tx signal c in one axis can be expressed by Equation 3 below.

$\begin{matrix} {{MathFigure}{\mspace{11mu} \;}3} & \; \\ {{{I(c)} = \frac{{- \left( {\sqrt{M_{n}} + 1 - {2i}} \right)}D_{n}}{2}},{i = 1},2,\ldots \mspace{14mu},\sqrt{M_{n}}} & \left\lbrack {{Math}.\mspace{14mu} 3} \right\rbrack \end{matrix}$

From the perspective of the I-axis, A(c) denotes a priori probability in an arbitrary Tx signal c. The Tx signal is c=(i_(n,1), i_(n,2), i_(n,3), . . . , i_(n,Bn)). In addition, i_(n,1) denotes a most significant bit (MSB), and i_(n,Bn) denotes a least significant bit (LSB). If y_(n) denotes a reception (Rx) signal, a log-likelihood ratio (LLR) of a bit i_(n,b) can be expressed by Equation 4 below.

$\begin{matrix} {{MathFigure}{\mspace{11mu} \;}4} & \; \\ {{{{LLR}\left( {i_{n,b}y_{n}} \right)} = {{\log \left\{ \frac{\sum\limits_{c \in {\{{i_{n,b} = 1}\}}}\; {{A(c)} \cdot {P\left( {y_{n}{I(c)}} \right)}}}{\sum\limits_{c \in {\{{i_{n,b} = 0}\}}}\; {{A(c)} \cdot {P\left( {y_{n}{I(c)}} \right)}}} \right\}} \approx \frac{\begin{pmatrix} {{2 \cdot y_{n} \cdot \left( {{I_{\min,1}(c)} - {I_{\min,0}(c)}} \right)} +} \\ \left( {{I_{\min,0}(c)}^{2} - {I_{\min,1}(c)}^{2} - {I_{\min,1}(c)}^{2}} \right) \end{pmatrix}}{\left( N_{0} \right)}}},{{{where}\mspace{14mu} {I_{\min,x}(c)}} = {\arg {\min\limits_{I{(c)}}\left( {y_{n} - {I\left( {c \in \left\{ {i_{n,b} = x} \right\}} \right)}} \right)}}}} & \left\lbrack {{Math}.\mspace{14mu} 4} \right\rbrack \end{matrix}$

An Rx signal z_(n) denotes a signal received from an n-th antenna and completely restored using a MIMO equalizer. If it is assumed that the Rx signal z creates an Rx signal y_(n) by removing channel information h_(n) through a zero-forcing (ZF) equalizer, the Rx signal y_(n) can be expressed as shown:

MathFigure 5

y _(n) =x _(n) +n′  [Math.5]

where x_(n) denotes a Tx signal. n′ denotes an increased noise generated when passing the ZF equalizer, and has a normal distribution of N(0, N₀/(2|h₀|²).

According to Equations 4 and 5, an average and dispersion of a received LLR value can be calculated as follows.

$\begin{matrix} {{MathFigure}{\mspace{11mu} \;}6} & \; \\ {{{{LLR}\left( {i_{n,b}y_{n}} \right)} \approx \frac{\begin{matrix} {{2 \cdot {E\left\lbrack x_{n} \right\rbrack} \cdot \left( {{I_{\min,1}(c)} - {I_{\min,0}(c)}} \right)} +} \\ \left( {{I_{\min,0}(c)}^{2} - {I_{\min,1}(c)}^{2}} \right) \end{matrix}}{N_{0}/{h_{n}}^{2}}}{{V\left\lbrack {{LLR}\left( {i_{n,b}y_{n}} \right)} \right\rbrack} \approx {\left( \frac{2 \cdot \left( {{I_{\min,1}(c)} - {I_{\min,0}(c)}} \right)}{N_{0}/{h_{n}}^{2}} \right)^{2} \cdot \frac{N_{0}}{2{h_{n}}^{2}}}}} & \left\lbrack {{Math}.\mspace{14mu} 6} \right\rbrack \end{matrix}$

If it is assumed that LLR values are calculated by performing MIMO combining with a bit unit in every retransmission, the LLR values are summed until rth retransmission is performed. The obtained sum of LLR values can be expressed by the following equation.

$\begin{matrix} {{MathFigure}{\mspace{11mu} \;}7} & \; \\ {{{{LLR}^{({comb})}\left( {i_{n,b}y_{n}} \right)} = {\sum\limits_{i = 0}^{r}\; {{LLR}^{(l)}\left( {i_{n,b}y_{n}} \right)}}},{b = 1},2,\ldots \mspace{14mu},B_{n}} & \left\lbrack {{Math}.\mspace{14mu} 7} \right\rbrack \end{matrix}$

By using a Q-function and a Chernoff bound, a bit error rate (BER) value of individual transmission bits constituting each data symbol can be expressed by the following equation.

$\begin{matrix} {{MathFigure}{\mspace{11mu} \;}8} & \; \\ \begin{matrix} {P_{n,{bc}}^{(r)} = {Q\left( \frac{{E\left\lbrack {{LLR}^{({comb})}\left( {i_{n,b}y_{n}} \right)} \right\rbrack}}{\sqrt{V\left\lbrack {{LLR}^{({comb})}\left( {i_{n,b}y_{n}} \right)} \right\rbrack}} \right)}} \\ {= {{Q\left( \frac{\sum\limits_{l = 0}^{r}\; {{E\left\lbrack {{LLR}^{(l)}\left( {i_{n,b}y_{n}} \right)} \right\rbrack}}}{\sqrt{\sum\limits_{l = 0}^{r}\; {{V\left\lbrack {{LLR}^{(l)}\left( {i_{n,b}y_{n}} \right)} \right\rbrack}}}} \right)} \leq}} \\ {{\frac{1}{2}\exp \left\{ {{- \frac{1}{2}} \cdot \frac{\left( {\sum\limits_{l = 0}^{r}\; {{E\left\lbrack {{LLR}^{(l)}\left( {i_{n,b}y_{n}} \right)} \right\rbrack}}} \right)^{2}}{\sum\limits_{l = 0}^{r}\; {V\left\lbrack {{LLR}^{(l)}\left( {i_{n,b}y_{n}} \right)} \right\rbrack}}} \right\}}} \\ {= {\frac{1}{2}\exp \left\{ {{- \frac{1}{2}} \cdot R_{n,{bc}}^{(r)}} \right\}}} \end{matrix} & \left\lbrack {{Math}.\mspace{14mu} 8} \right\rbrack \end{matrix}$

Herein, R^((r)) _(n,b|c) denotes a ratio of an average and a dispersion of the summed LLR values expressed in a last part of Equation 8, and is defined as a conditional bit reliability of a bit i_(n,b). The term “conditional” means that the bit reliability is at a time when a Tx signal (i.e., data symbol) c is transmitted. By averaging conditional bit reliabilities for all possible data symbols, an average bit reliability can be obtained by the following equation.

$\begin{matrix} {{MathFigure}{\mspace{11mu} \;}9} & \; \\ {R_{n,b}^{(r)} = {\frac{1}{\sqrt{M_{n}}}{\sum\limits_{c}\; R_{n,{bc}}^{(r)}}}} & \left\lbrack {{Math}.\mspace{14mu} 9} \right\rbrack \end{matrix}$

According to Equation 9, the conditional bit reliability and the average bit reliability depending on a transmission bit of the M-QAM can be obtained.

Table 1 shows the average bit reliability in the I-axis of 16-QAM when initial transmission is performed. Herein, D² _(n)=E² _(s)/10.

TABLE 1 c I(c) R_(n, 1c)^((r)) R_(n, 2c)^((r)) 00 −3D_(n)/2 $2 \cdot \frac{D_{n}^{2}{h_{n}}^{2}}{N_{0}}$ $\frac{1}{2} \cdot \frac{D_{n}^{2}{h_{n}}^{2}}{N_{0}}$ 01 −D_(n)/2 $\frac{1}{2} \cdot \frac{D_{n}^{2}{h_{n}}^{2}}{N_{0}}$ $\frac{1}{2} \cdot \frac{D_{n}^{2}{h_{n}}^{2}}{N_{0}}$ 10 3D_(n)/2 $2 \cdot \frac{D_{n}^{2}{h_{n}}^{2}}{N_{0}}$ $\frac{1}{2} \cdot \frac{D_{n}^{2}{h_{n}}^{2}}{N_{0}}$ 11 D_(n)/2 $\frac{1}{2} \cdot \frac{D_{n}^{2}{h_{n}}^{2}}{N_{0}}$ $\frac{1}{2} \cdot \frac{D_{n}^{2}{h_{n}}^{2}}{N_{0}}$ R_(n, b)^((r)) N/A $1.25 \cdot \frac{D_{n}^{2}{h_{n}}^{2}}{N_{0}}$ $\frac{1}{2} \cdot \frac{D_{n}^{2}{h_{n}}^{2}}{N_{0}}$

Table 2 shows the average bit reliability in the I-axis of 64-QAM when initial transmission is performed. Herein, D² _(n)=E² _(s)/84.

TABLE 2 c I(c) R_(n, 1c)^((r)) R_(n, 2c)^((r)) R_(n, 3c)^((r)) 000 −7D_(n)/2 $8 \cdot \frac{D_{n}^{2}{h_{n}}^{2}}{N_{0}}$ $2 \cdot \frac{D_{n}^{2}{h_{n}}^{2}}{N_{0}}$ $\frac{1}{2} \cdot \frac{D_{n}^{2}{h_{n}}^{2}}{N_{0}}$ 001 −5D_(n)/2 $4.5 \cdot \frac{D_{n}^{2}{h_{n}}^{2}}{N_{0}}$ $\frac{1}{2} \cdot \frac{D_{n}^{2}{h_{n}}^{2}}{N_{0}}$ $\frac{1}{2} \cdot \frac{D_{n}^{2}{h_{n}}^{2}}{N_{0}}$ 010 −D_(n)/2 $\frac{1}{2} \cdot \frac{D_{n}^{2}{h_{n}}^{2}}{N_{0}}$ $2 \cdot \frac{D_{n}^{2}{h_{n}}^{2}}{N_{0}}$ $\frac{1}{2} \cdot \frac{D_{n}^{2}{h_{n}}^{2}}{N_{0}}$ 011 −3D_(n)/2 $2 \cdot \frac{D_{n}^{2}{h_{n}}^{2}}{N_{0}}$ $\frac{1}{2} \cdot \frac{D_{n}^{2}{h_{n}}^{2}}{N_{0}}$ $\frac{1}{2} \cdot \frac{D_{n}^{2}{h_{n}}^{2}}{N_{0}}$ 100 7D_(n)/2 $8 \cdot \frac{D_{n}^{2}{h_{n}}^{2}}{N_{0}}$ $2 \cdot \frac{D_{n}^{2}{h_{n}}^{2}}{N_{0}}$ $\frac{1}{2} \cdot \frac{D_{n}^{2}{h_{n}}^{2}}{N_{0}}$ 101 5D_(n)/2 $4.5 \cdot \frac{D_{n}^{2}{h_{n}}^{2}}{N_{0}}$ $\frac{1}{2} \cdot \frac{D_{n}^{2}{h_{n}}^{2}}{N_{0}}$ $\frac{1}{2} \cdot \frac{D_{n}^{2}{h_{n}}^{2}}{N_{0}}$ 110 D_(n)/2 $\frac{1}{2} \cdot \frac{D_{n}^{2}{h_{n}}^{2}}{N_{0}}$ $2 \cdot \frac{D_{n}^{2}{h_{n}}^{2}}{N_{0}}$ $\frac{1}{2} \cdot \frac{D_{n}^{2}{h_{n}}^{2}}{N_{0}}$ 111 3D_(n)/2 $2 \cdot \frac{D_{n}^{2}{h_{n}}^{2}}{N_{0}}$ $\frac{1}{2} \cdot \frac{D_{n}^{2}{h_{n}}^{2}}{N_{0}}$ $\frac{1}{2} \cdot \frac{D_{n}^{2}{h_{n}}^{2}}{N_{0}}$ R_(n, b)^((r)) N/A $3.75 \cdot \frac{D_{n}^{2}{h_{n}}^{2}}{N_{0}}$ $1.25 \cdot \frac{D_{n}^{2}{h_{n}}^{2}}{N_{0}}$ $\frac{1}{2} \cdot \frac{D_{n}^{2}{h_{n}}^{2}}{N_{0}}$

As shown in Tables 1 and 2, the average bit reliability is higher in the MSB (i.e., i_(n,1)) than in the LSB (i.e., i_(n,2)). This implies that more bit errors can be resulted when the same information is transmitted through an LSB position than when the same information is transmitted through an MSB position. Therefore, if a bit transmitted through the LSB in initial transmission can be transmitted through the MSB in retransmission, reliability of an overall link performance can be increased.

Comparing Tables 1 and 2, the average bit reliability has a significant difference according to a modulation scheme. Therefore, the link performance can be more improved when bit swapping is performed between data symbols using different modulation schemes than when bit swapping is performed between data symbols using one modulation scheme in every retransmission. Swapping is defined as positional changes between bits.

Since all antennas use the same modulation scheme in single-modulation transmission, data symbols for all transmit antennas have the same type of average bit reliability. Therefore, theoretically, the average bit reliability of all data symbols can be determined to be almost constant by performing several retransmissions. However, since data symbols for all transmit antennas have different types of average bit reliability in multi-modulation transmission, swapping has to be performed so that a difference of average bit reliability can be reduced as much as possible.

To quantize the difference of average bit reliability by performing bit swapping, a minimum squared error (MSE) is defined by the following equation.

$\begin{matrix} {{MathFigure}{\mspace{11mu} \;}10} & \; \\ {ɛ^{(r)} = {\sum\limits_{n = 1}^{N_{t}}\; {\sum\limits_{b = 1}^{B_{n}}\; \left( {R_{n,b}^{(r)} - {E\left\lbrack R_{n,b}^{(r)} \right\rbrack}} \right)^{2}}}} & \left\lbrack {{Math}.\mspace{14mu} 10} \right\rbrack \end{matrix}$

By minimizing the MSE of bit reliability defined in Equation 10, the difference of average reliability can be reduced in every retransmission. Each bit has the difference of average reliability in multi-modulation transmission.

In R^((r−1)) _(n,b) obtained by performing (r−1) retransmissions, a best swapping matrix I(r) is a matrix that satisfies the following equation.

$\begin{matrix} {{MathFigure}{\mspace{11mu} \;}11} & \; \\ {I^{(r)} = {\arg {\min\limits_{\forall I^{(r)}}ɛ^{(r)}}}} & \left\lbrack {{Math}.\mspace{14mu} 11} \right\rbrack \end{matrix}$

If the number of transmit antennas and a modulation scheme used by each transmit antenna are predetermined, an average bit reliability to be added to an LLR value for each bit in every retransmission is predetermined. Therefore, when an average bit reliability for each bit is determined after performing (r−1) retransmissions, a bit having a lowest bit reliability is shifted to a position at which a highest bit reliability can be obtained in rth retransmission, and a bit having a highest bit reliability is shifted to a position at which a lowest bit reliability can be obtained in rth retransmission. Remaining bits are also shifted to positions having a high reliability in an ascending order of bit reliability. According to this method, conditions for minimizing the MSE of Equation 10 is satisfied in a state that an average bit reliability value summed in every transmission is determined for each transmit antenna and for each bit.

A minimum size of a bit swapping set according to the proposed method is defined as a first η satisfying ε^((η−1)<ε) ^((η).)

For example, it is assumed that a 1st data symbol modulated with 16-QAM is transmitted through a 1st transmit antenna, and a 2nd data symbol modulated with 64-QAM is transmitted through a 2nd transmit antenna.

FIG. 6 is a graph showing changes in an MSE of average bit reliability with respect to the number of retransmissions according to the proposed bit swapping scheme. The graph shows changes of ε^((r)) that is normalized to ε⁽⁰⁾ of initial transmission.

Referring to FIG. 6, according to the proposed bit swapping scheme, a difference of bit reliability can be reduced by half in 1st retransmission. In this case, a minimum size of the bit swapping set is 3. Therefore, the bit swapping set proposed in this example can be expressed by the following equation.

$\begin{matrix} {{MathFigure}{\mspace{11mu} \;}12} & \; \\ {I^{(r)} \in \begin{Bmatrix} {\begin{bmatrix} i_{1,1} & i_{1,2} & i_{1,3} \\ i_{2,1} & i_{2,2} & x \end{bmatrix},} \\ {\begin{bmatrix} i_{1,1} & i_{2,2} & i_{2,1} \\ i_{1,3} & i_{i,2} & x \end{bmatrix},} \\ {\begin{bmatrix} i_{1,1} & i_{1,3} & i_{2,1} \\ i_{1,2} & i_{2,2} & x \end{bmatrix},} \end{Bmatrix}} & \left\lbrack {{Math}.\mspace{14mu} 12} \right\rbrack \end{matrix}$

In each matrix, a first row represents 64-QAM transmission, and a second row represents 16-QAM transmission. A first element of the swapping set represents initial transmission. As shown in Tables 1 and 2 above, in case of 64-QAM and 16-QAM, an MSB of 16-QAM is positioned where a highest bit reliability can be obtained, followed by an LSB of 16-QAM, a middle bit of 64-QAM, and an LSB of 64-QAM, in that order. Therefore, a second matrix is configured such that an LSB of 64-QAM which has a lowest bit reliability in a first matrix is positioned to an MSB of 16-QAM, and a middle bit of 64-QAM is positioned to an LSB of 16-QAM. Accordingly, an MSE of average bit reliability is reduced. In a third matrix, bits are mapped by reverse-sorting according to the summed average bit reliability obtained by performing two retransmissions. Since the average bit reliability is sufficiently decreased after performing three retransmissions as shown in FIG. 6, the proposed bit mapping of Equation 12 is repeatedly used.

A method of reducing a difference of average bit reliability between bits by performing bit swapping is necessary to find a best BSI scheme. However, it is difficult to reduce a difference of conditional bit reliability between data symbols when only bit swapping is used.

To reduce the difference of conditional bit reliability, it is necessary to perform an inversion operation in retransmission. In a case where (1,0) and (1,1) of Table 1 are transmitted, if an LSB is inversely mapped in transmission, a receiver receives (1,1) and (1,0) having an opposite LLR value. When using single modulation transmission, since an MSB of 16-QAM has a total of two absolute values for different LLR values, an absolute value difference of LLR between data symbols can be reduced when transmission is made by performing the inversion operation one time for initial transmission. However, when using multi-modulation transmission, rows of a Tx block are configured to have different sizes. Bit inversion for an nth antenna whose row is less than B_(max) has to be repetitively performed. Bit inversion on an n-th row of the Tx block can be expressed by the following equation.

$\begin{matrix} {{MathFigure}{\mspace{11mu} \;}13} & \; \\ {{I^{(r)}\left( {n,\text{:}} \right)} \in \begin{Bmatrix} \begin{bmatrix} i_{n,1} & \overset{\_}{i_{n,2}} & i_{n,3} & \ldots & i_{n,B_{n}} \end{bmatrix} \\ \begin{bmatrix} i_{n,1} & i_{n,2} & \overset{\_}{i_{n,3}} & \ldots & i_{n,B_{n}} \end{bmatrix} \\ \begin{matrix} \; & \; & \vdots & \; & \; \end{matrix} \\ \begin{bmatrix} i_{n,1} & \overset{\_}{i_{n,2}} & \overset{\_}{i_{n,3}} & \ldots & \overset{\_}{i_{n,B_{n}}} \end{bmatrix} \\ \begin{matrix} \; & \; & \vdots & \; & \; \end{matrix} \\ \begin{bmatrix} i_{n,1} & \overset{\_}{i_{n,2}} & i_{n,3} & \ldots & i_{n,B_{n}} \end{bmatrix} \\ \begin{bmatrix} i_{n,1} & i_{n,2} & \overset{\_}{i_{n,3}} & \ldots & i_{n,B_{n}} \end{bmatrix} \\ \begin{matrix} \; & \; & \vdots & \; & \; \end{matrix} \\ \begin{bmatrix} i_{n,1} & \overset{\_}{i_{n,2}} & \overset{\_}{i_{n,3}} & \ldots & \overset{\_}{i_{n,B_{n}}} \end{bmatrix} \end{Bmatrix}} & \left\lbrack {{Math}.\mspace{14mu} 13} \right\rbrack \end{matrix}$

In the above equation, I^((r))(n,:) denotes an nth row of I^((r)). A size of the proposed inversion matrix is 2^(Bmax−1). The number of times of repeating the inversion operation within I^((r))(n,:) is 2^(Bmax−Bn). According to the principal of the above inversion operation, the bit reliability exactly the same for all data symbols in multi-modulation transmission.

By combining bit swapping and bit inversion, combined bit reliabilities can be approximated with respect to all data symbols and all bit positions. Therefore, a BSI set size proposed in the present invention for multi-modulation transmission is η2^(Bmax−1).

FIG. 7 shows a BSI scheme in a 2×2 MIMO system using 16-QAM and 64-QAM according to an embodiment of the present invention. Herein, Equations 12 and 13 are used in combination by the 2×2 MIMO system using 16-QAM and 64-QAM.

Referring to FIG. 7, three Tx blocks are obtained by performing bit swapping according to Equation 12 above. If the three Tx blocks are extended by performing an inversion operation, a total of 12 BSI sets are obtained. A total of four Tx blocks are configured by performing the inversion operation with respect to a 1st transmit antenna using 64-QAM. Four Tx blocks are configured by repeatedly performing an inversion and non-inversion operation two times on a second bit with respect to a 2nd transmit antenna using 16-QAM whose modulation order is lower than that of 64-QAM. A difference of average bit reliability can be significantly reduced using the inversion operation.

Now, optimization of the BSA scheme in multi-modulation transmission will be described.

In a diversity transmission scheme used in a conventional MIMO system, the same data symbol is transmitted to a plurality of transmit antennas to obtain spatial diversity, and the same data symbol is repeatedly transmitted at different time to obtain time diversity. In the HARQ scheme, a time difference between retransmissions is several-millisecond (ms) unit, which is significantly greater than a symbol unit. Thus, a channel condition is significantly different between retransmissions. Accordingly, when retransmission is considered, a higher diversity gain can be obtained using the spatial diversity in which a transmit antenna is changed for each bit.

FIG. 8 shows an example of a BSA scheme in single modulation transmission.

Referring to FIG. 8, transmission is achieved by applying different shift amounts to respective columns so that bits constituting one data symbol are transmitted respectively through different antennas in every retransmission. Applying different shift amounts to respective columns of a Tx block is called multi-step cyclic shift (MSCS). A maximum diversity gain can be expected through the MSCS.

Since respective bits constituting one data symbol can be transmitted through different transmit antennas, if the number B of bits constituting one I-axis or Q-axis symbol is greater than the number Nt of transmit antennas, up to B additional diversities can be obtained.

The MSCS can be expressed with the aforementioned matrix expression according to the following equation.

MathFigure 14

I ⁽¹⁾ =f _(MSCS)(I ⁽⁰⁾)  [Math.14]

In case of performing two or more retransmissions, a time difference between retransmissions is relatively large, and thus it can be regarded that a channel response is independent between retransmissions. Therefore, it can be regarded that immediately previous BSA does not have an effect in determining of BSA in every retransmission. Accordingly, BSA for 2nd or higher retransmission as shown Equation 15 can also directly use an optimal BSA scheme of Equation 14 irrespective of previous BSA.

MathFigure 15

I ^((r)) =f _(MSCS)(I ^(r−1))), r=1, 2, . . . , r_(max)  [Math.15]

In multi-modulation transmission, the BSA in the aforementioned single modulation transmission can be extended in the similar manner. It is desired to use MSCS-based BSA since the purpose of using the multi-modulation transmission is to obtain a high spatial diversity by transmitting bits constituting one data symbol through different transmit antennas. However, since the number of bits constituting each row of a Tx block is different unlike in the single modulation transmission, a group MSCS (GMSCS) scheme is proposed.

According to the GMSCS scheme, vertical swapping of bits through MSCS is performed only between rows having the same number of bits. If a given matrix is divided into an inter-row permutation matrix P_(row) and a depermutation matrix D_(row), an operation of GMSCS can be mathematically expressed as follows.

$\begin{matrix} {{MathFigure}{\mspace{11mu} \;}16} & \; \\ \begin{matrix} {I^{(r)} = {{f_{GMSCS}\left( {I^{({r - 1})} \cdot P_{row}} \right)} \cdot D_{row}}} \\ {= {{f_{GMSCS}\left( \begin{bmatrix} I_{2}^{({r - 1})} \\ I_{3}^{({r - 1})} \\ \vdots \\ I_{B_{\max}}^{({r - 1})} \end{bmatrix} \right)} \cdot D_{row}}} \\ {{= {\begin{bmatrix} {f_{MSCS}\left( I_{2}^{({r - 1})} \right)} \\ {f_{MSCS}\left( I_{3}^{({r - 1})} \right)} \\ \vdots \\ {f_{MSCS}\left( I_{B_{\max}}^{({r - 1})} \right)} \end{bmatrix} \cdot D_{row}}},{r = 1},2,\ldots \mspace{14mu},r_{\max}} \end{matrix} & \left\lbrack {{Math}.\mspace{14mu} 16} \right\rbrack \end{matrix}$

Herein, I^((r−1)) _(m) denotes a matrix which is newly created by collecting rows having m bits. The operation of BSA can maintain independency from the operation of BSI. That is, an LLR value equally adjusted by the BSI scheme is not changed.

FIG. 9 shows a BSA scheme in a 5×5 MIMO system using 16-QAM and 64-QAM according to an embodiment of the present invention.

Referring to FIGS. 9, 1st, 3rd, and 5th transmit antennas use 64-QAM, and 2nd and 4th antennas use 16-QAM. Therefore, in initial transmission, a Tx block is expressed in a format of I⁽⁰⁾.

A matrix I⁽¹⁾ is generated by an algorithm using GMSCS. When bit swapping is performed between 2nd and 4th transmit antennas using 16-QAM, a second bit of each row is swapped by 1-step cyclic shift. When bit swapping is performed between 1st, 3rd, and 5th transmit antennas using 64-QAM, a second bit is swapped by 1-step cyclic shift, and a third bit is swapped by 2-step cyclic shift. Matrixes I⁽²⁾, . . . , I^((r)) for two or more retransmissions can also be obtained using the same GMSCS algorithm.

FIG. 10 shows an example of applying a BSI scheme and a BSA scheme to a 2×2 MIMO system using 16-QAM and 64-QAM. A combination of the BSI and BSA schemes is used in multi-modulation transmission. In this example, the matrix set and the BSI scheme of FIG. 7 are applied. This is because while the proposed BSA scheme using GMSCS is applied between two or more transmit antennas using the same modulation scheme, the proposed BSA scheme is not applied in this example since there are only two transmit antennas using different modulation schemes.

The proposed scheme is applied for the number of arbitrary retransmissions, a combination of arbitrary m-QAM schemes, and the number of arbitrary transmit antennas.

FIG. 11 is a graph showing a result obtained by performing simulations in a 2×2 MIMO system using 16-QAM and 64-QAM. The result of FIG. 11 is obtained by performing one retransmission with respect to the conventional method, an antenna switching method in which only BSA is performed in a symbol unit, and the proposed method. It is shown that a gain of about 3 dB can be obtained in the proposed method in comparison with the conventional method irrespective of a speed of a UE.

II. Constellation Rearrangement for MIMO Systems with HARQ

Spatial Multiplexing (SM) system using N transmit antenna and M receive antennas is considered. Suppose that the information bits are encoded with a coding rate of R by turbo code and L-coded bits are mapped into one of the complex constellation points out of the 2L-ary Quadrature Amplitude Modulation (QAM) symbols. The turbo code is used as exemplary purpose only, any well-known coding scheme such as convolution code, etc. may be used. Each symbol is modulated by Inverse Fast Fourier Transform (IFFT). In the transmitter, every N OFDM symbols can be transmitted through N transmit antennas at the same time.

From now on, we will analyze for 2L-QAM symbol in terms of unequal error probability over all of component bits. As known well, in square constellation such as 2^(L)-QAM symbol, there exist symmetries between In-phase and Quadrature components. Therefore, in order to make analysis simpler, we just consider one dimension 4-PAM constellation for either In-phase or Quadrature component of 16-QAM.

FIG. 12 shows Gray mapping for 16-QAM and FIG. 13 shows decision boundary to calculate the probability of bit error. It is assumed that 16-QAM is modulated by gray mapping and 4-PAM constellation for Most Significant Bit (MSB) and Least Significant Bit (LSB) to calculate the probability of bit error on each component bit is subject to one dimension, either In-phase or Quadrature. All possible 4-PAM symbols are s₁, s₂, s₃ and s₄. Also these symbols are located at −d₁, −d₂, −d₁, +d₂, d₁, −d₂ and d₁, +d₂, respectively, in which 2d₁ means the distance between two decision boundary within 4-PAM constellation for LSB and 2d₂ represents the distance between two neighboring symbols within the 4-PAM constellation.

The probability of bit error for MSB, which is the probability of bit error when s1 symbol (i.e. 10) is transmitted, can be expressed as shown:

$\begin{matrix} {{MathFigure}\mspace{14mu} 17} & \; \\ \begin{matrix} {P_{e} = {P\left\{ {{{\langle{n,\phi}\rangle} \geq \frac{d}{2}} = \frac{2\left( {d_{1} + d_{2}} \right)}{2}} \right\}}} \\ {= {\int_{d_{1} + d_{2}}^{\infty}{\frac{1}{\sqrt{2{\pi\sigma}_{n}^{2}}}^{\frac{1}{2\sigma_{n}^{2}}u^{2}}\ {u}}}} \\ {= {\int_{\frac{d_{1} + d_{2}}{\sigma_{n}^{2}}}^{\infty}{\frac{1}{\sqrt{2\pi}}^{\frac{1}{2}u^{2}}\ {u}}}} \\ {= {Q\left( \frac{d_{1} + d_{2}}{\sigma_{n}^{2}} \right)}} \end{matrix} & \left\lbrack {{Math}.\mspace{14mu} 17} \right\rbrack \end{matrix}$

where n means white Gaussian noise with variance σ² _(n), and Ψ is a unit norm vector along the line between s₁ and s₄. <x, y> represents the inner product between x and y. The probability of bit error can be defined as the probability that the component of the noise vector n along the line connecting between the two associated symbols is greater than half the distance along this line. The error probability of s₄ is the same as that of s₁. In the same way, when s₂ and s₃ is transmitted, the probability of bit error is given as shown below.

$\begin{matrix} {{MathFigure}\mspace{14mu} 18} & \; \\ {Q\left( \frac{d_{1} + d_{2}}{\sigma_{n}} \right)} & \left\lbrack {{Math}.\mspace{14mu} 18} \right\rbrack \end{matrix}$

Thus, the average probability of bit error for MSB can be obtained as shown:

$\begin{matrix} {{MathFigure}\mspace{14mu} 19} & \; \\ \begin{matrix} {P_{e} = {{P_{e{s\; 1}}{p\left( s_{1} \right)}} + {P_{e{s\; 2}}{p\left( s_{2} \right)}} + {P_{e{s\; 3}}{p\left( s_{3} \right)}} + {P_{e{s\; 4}}{p\left( s_{4} \right)}}}} \\ {= {{\frac{1}{4}{Q\left( \frac{d_{1} + d_{2}}{\sigma_{n}} \right)}} + {\frac{1}{4}{Q\left( \frac{d_{1} - d_{2}}{\sigma_{n}} \right)}} + {\frac{1}{4}{Q\left( \frac{d_{1} - d_{2}}{\sigma_{n}} \right)}} +}} \\ {{\frac{1}{4}{Q\left( \frac{d_{1} + d_{2}}{\sigma_{n}} \right)}}} \end{matrix} & \left\lbrack {{Math}.\mspace{14mu} 19} \right\rbrack \end{matrix}$

where p(s_(i)) represents the priori probability of s_(i) and P_(e|si) is the probability of bit error when s_(i) is transmitted to the receiver. For the second, the probability for LSB as well as the way to obtain the probability of bit error for MSB can be calculated. The probability of bit error for LSB is thereby given as shown below.

$\begin{matrix} {{MathFigure}\mspace{14mu} 20} & \; \\ {{\frac{1}{4}\left( {{Q\left( \frac{d_{2}}{\sigma_{n}} \right)} - {Q\left( \frac{{2d_{1}} + d_{2}}{\sigma_{n}} \right)}} \right)} + {\frac{1}{4}\left( {{Q\left( \frac{d_{2}}{\sigma_{n}} \right)} + {Q\left( \frac{{2d_{1}} + d_{2}}{\sigma_{n}} \right)}} \right)} + {\frac{1}{4}\left( {{Q\left( \frac{d_{2}}{\sigma_{n}} \right)} + {Q\left( \frac{{2d_{1}} + d_{2}}{\sigma_{n}} \right)}} \right)} + {\frac{1}{4}\left( {{Q\left( \frac{d_{2}}{\sigma_{n}} \right)} - {Q\left( \frac{{2d_{1}} + d_{2}}{\sigma_{n}} \right)}} \right)}} & \left\lbrack {{Math}.\mspace{14mu} 20} \right\rbrack \end{matrix}$

Now, it can be known from Equations 19 and 20 that the component bits in either Inphase or Quadrature have the different probability of bit error. The difference of probability of bit error between all possible candidate bits to be either MSB or LSB can be also observed. Furthermore, when the probability of bit error is calculated, because there are two decision boundaries for LSB, the decision boundary is easy to be influenced by the wrong information from the other side as shown in Equation 20.

Also, let dd₁=2d₂=2d, we can approximate Equation 19 as follows.

$\begin{matrix} {{MathFigure}\mspace{14mu} 21} & \; \\ {{\frac{1}{4}{Q\left( \frac{3d}{\sigma_{n}} \right)}} + {\frac{1}{4}{Q\left( \frac{d}{\sigma_{n}} \right)}} + {\frac{1}{4}{Q\left( \frac{d}{\sigma_{n}} \right)}} + {\frac{1}{4}{Q\left( \frac{3d}{\sigma_{n}} \right)}}} & \left\lbrack {{Math}.\mspace{14mu} 21} \right\rbrack \end{matrix}$

Similarly, assuming that the effect from the other decision boundary can be ignored, the probability of bit error for LSB can be represented as shown below.

$\begin{matrix} {{MathFigure}\mspace{14mu} 22} & \; \\ {{\frac{1}{4}{Q\left( \frac{d}{\sigma_{n}} \right)}} + {\frac{1}{4}{Q\left( \frac{d}{\sigma_{n}} \right)}} + {\frac{1}{4}{Q\left( \frac{d}{\sigma_{n}} \right)}} + {\frac{1}{4}{Q\left( \frac{d}{\sigma_{n}} \right)}}} & \left\lbrack {{Math}.\mspace{14mu} 22} \right\rbrack \end{matrix}$

For HARQ-Chase Combing (CC), it is assumed that the probability of bit error for the case that the current Log Likelihood Ratio (LLR) values and the previous LLR values are combined in a demapper. According to the type of channel environment, the probability of error for HARQ-CC into three sorts can be divided as follows.

a) Additive White Gaussian Noise (AWGN) channel:

$\begin{matrix} {{MathFigure}\mspace{14mu} 23} & \; \\ {Q\left( \sqrt{\frac{1}{2\sigma_{n}^{2}}{\sum\limits_{i = 1}^{M}\; {d_{i}}^{2}}} \right)} & \left\lbrack {{Math}.\mspace{14mu} 23} \right\rbrack \end{matrix}$

where d denotes a Euclidean distance at i-th retransmission, and M means the maximum number of retransmissions.

b) Fading channel with incoherent demodulation:

$\begin{matrix} {{MathFigure}\mspace{14mu} 24} & \; \\ {Q\left( \sqrt{\frac{1}{2\sigma_{n}^{2}}{\sum\limits_{i = 1}^{M}\; {{h_{i}}^{2}{d_{i}}^{2}}}} \right)} & \left\lbrack {{Math}.\mspace{14mu} 24} \right\rbrack \end{matrix}$

where h means the channel value at i-th retransmission. For a receiver with the incoherent demodulation, the probability of error of HARQ-CC in fading channel environment is given by Equation 24.

c) Fading channel with coherent demodulation:

$\begin{matrix} {{MathFigure}\mspace{14mu} 25} & \; \\ {Q\left( \sqrt{\frac{1}{2\sigma_{n}^{2}}{\sum\limits_{i = 1}^{M}\; \left\lbrack {{{{{Re}\left\{ h_{i} \right\}}}^{2}{{{Re}\left\{ d_{i} \right\}}}^{2}} + {{{{Im}\left\{ h_{i} \right\}}}^{2}{{{Im}\left\{ d_{i} \right\}}}^{2}}} \right\rbrack}} \right)} & \left\lbrack {{Math}.\mspace{14mu} 25} \right\rbrack \end{matrix}$

The probability of error for HARQ-CC in fading channel environment with coherent demodulation is given by Equation 25.

Therefore, we consider both the probability of error for HARQ-CC and the probability of bit error for MSB and LSB in the course of chase combining as shown in Table 3.

TABLE 3 component Initial Transmission Retransmission (CC) bits MSB LSB MSB LSB 10 $\frac{1}{4}{Q\left( \frac{3d}{\sigma_{n}^{2}} \right)}$ $\frac{1}{4}{Q\left( \frac{d}{\sigma_{n}^{2}} \right)}$ $\frac{1}{4}{Q\left( \frac{3\sqrt{2}d}{\sigma_{n}} \right)}$ $\frac{1}{4}{Q\left( \frac{\sqrt{2}d}{\sigma_{n}^{2}} \right)}$ 11 $\frac{1}{4}{Q\left( \frac{d}{\sigma_{n}^{2}} \right)}$ $\frac{1}{4}{Q\left( \frac{d}{\sigma_{n}^{2}} \right)}$ $\frac{1}{4}{Q\left( \frac{\sqrt{2}d}{\sigma_{n}} \right)}$ $\frac{1}{4}{Q\left( \frac{\sqrt{2}d}{\sigma_{n}^{2}} \right)}$ 01 $\frac{1}{4}{Q\left( \frac{d}{\sigma_{n}^{2}} \right)}$ $\frac{1}{4}{Q\left( \frac{d}{\sigma_{n}^{2}} \right)}$ $\frac{1}{4}{Q\left( \frac{\sqrt{2}d}{\sigma_{n}} \right)}$ $\frac{1}{4}{Q\left( \frac{\sqrt{2}d}{\sigma_{n}^{2}} \right)}$ 00 $\frac{1}{4}{Q\left( \frac{3d}{\sigma_{n}^{2}} \right)}$ $\frac{1}{4}{Q\left( \frac{d}{\sigma_{n}^{2}} \right)}$ $\frac{1}{4}{Q\left( \frac{3\sqrt{2}d}{\sigma_{n}} \right)}$ $\frac{1}{4}{Q\left( \frac{\sqrt{2}d}{\sigma_{n}^{2}} \right)}$ Average ${\frac{1}{2}{Q\left( \frac{3d}{\sigma_{n}^{2}} \right)}} + {\frac{1}{2}{Q\left( \frac{d}{\sigma_{n}^{2}} \right)}}$ $Q\left( \frac{d}{\sigma_{n}^{2}} \right)$ ${\frac{1}{2}{Q\left( \frac{3\sqrt{2}d}{\sigma_{n}} \right)}} + {\frac{1}{2}{Q\left( \frac{\sqrt{2}d}{\sigma_{n}} \right)}}$ $Q\left( \frac{\sqrt{2}d}{\sigma_{n}^{2}} \right)$

Table 3 shows that the difference of bit error probability still exists during retransmission. Moreover, although a little improvement of Signal-to-Noise Ratio (SNR) can be obtained, the difference of bit error probability between MSB and LSB is further increased comparing to the case of initial transmission. In case of MSB, the difference of bit error probability between component bits which are subjected to MSB is also increased.

Horizontal bitwise optimization makes a role to reduce the difference of reliability between coded bits. The performance of decoding depends on how much the distribution of reliability of bits can be uniform. As the distribution of bits becomes to be uniform, more coding gain can be obtained. Therefore, horizontal bitwise operations are used to maximize coding gain in retransmission. We regard two operations such as swapping and inversion as horizontal bitwise mapping. The swapping operation makes the decision boundary corresponding to MSB and LSB to be exchanged each other during retransmission.

FIG. 14 shows swapping operation and FIG. 15 shows inversion operation. In FIG. 13, it is shown that the number of decision boundary for MSB and LSB is one and two, respectively. When employing swapping operation during retransmission, the decision boundary for MSB and LSB is swapped to result in decreasing the difference of the bit error probability between MSB and LSB as shown in FIG. 14. Inversion bitwise operation makes the difference of bit error probability of candidate component bits to be decreased. In this case, decision boundary is not exchanged between MSB and LSB. As shown in FIG. 15, inversion for LSB is only valid. Therefore, when we need to average the probability of bit error between candidate component bits subjected to LSB, i.e. for the applying the inversion operation to a component bit which has not more than two decision boundary, we need swapping operation with inversion operation.

Under employing HARQ, the proposed design criterion to optimize horizontal bitwise mapping for 2^(L)-QAM is represented as shown:

$\begin{matrix} {{MathFigure}\mspace{14mu} 26} & \; \\ {\left\{ \mu_{m} \right\} = {\arg \underset{\mu_{i} \in \mu}{Max}{\begin{matrix} {{\frac{1}{2^{{L/2} - 1}L}{\sum\limits_{k}^{2^{L/2}}\; {\sum\limits_{j}^{L/2}\; {Q\left( \sqrt{\frac{1}{2\sigma_{n}^{2}}{\sum\limits_{i = 1}^{m}\; {{\mu_{i}\left( d_{j,k} \right)}}^{2}}} \right)}}}} -} \\ {\frac{1}{2^{{L/2} - 1}L}{\sum\limits_{k}^{2^{L/2}}\; {\sum\limits_{j}^{L/2}\; {Q\left( \sqrt{\frac{1}{2\sigma_{n}^{2}}{\sum\limits_{i = 1}^{m - 1}\; {{\mu_{i}\left( d_{j,k} \right)}}^{2}}} \right)}}}} \end{matrix}}}} & \left\lbrack {{Math}.\mspace{14mu} 26} \right\rbrack \end{matrix}$

where d_(j,k) represents the distance of k-th candidate component bit subject to j-th component bit (i.e. MSB or LSB etc.). μ_(i) denotes a horizontal mapping at i-th retransmission. μ indicates the set of horizontal mapping. The number of combination that we can compose from bitwise swapping and inversion is _(L/2)C₂·(2^(L/2)−1). But, we can reduce the number of combination through the assumptions that bit inversion operation is only valid when there are more than two decision boundaries and we can know the difference of bit error probability between component bits from the approach-based Q-function. In case of 16-QAM, optimal horizontal mapping set can be calculated by comparisons of between 4 cases at every retransmission.

Table 4 shows a comparison of bit error probability of the swapping operation during first retransmission and that of conventional chase combining.

TABLE 4 component Chase combining (none + none) Chase combining (none + swapping) bits MSB LSB MSB LSB 10 $\frac{1}{4}{Q\left( \frac{3\sqrt{2}d}{\sigma_{n}} \right)}$ $\frac{1}{4}{Q\left( \frac{\sqrt{2}d}{\sigma_{n}^{2}} \right)}$ $\frac{1}{4}{Q\left( \frac{\sqrt{10}d}{\sigma_{n}} \right)}$ $\frac{1}{4}{Q\left( \frac{\sqrt{10}d}{\sigma_{n}} \right)}$ 11 $\frac{1}{4}{Q\left( \frac{\sqrt{2}d}{\sigma_{n}} \right)}$ $\frac{1}{4}{Q\left( \frac{\sqrt{2}d}{\sigma_{n}^{2}} \right)}$ $\frac{1}{4}{Q\left( \frac{\sqrt{2}d}{\sigma_{n}} \right)}$ $\frac{1}{4}{Q\left( \frac{\sqrt{2}d}{\sigma_{n}} \right)}$ 01 $\frac{1}{4}{Q\left( \frac{\sqrt{2}d}{\sigma_{n}} \right)}$ $\frac{1}{4}{Q\left( \frac{\sqrt{2}d}{\sigma_{n}^{2}} \right)}$ $\frac{1}{4}{Q\left( \frac{\sqrt{2}d}{\sigma_{n}} \right)}$ $\frac{1}{4}{Q\left( \frac{\sqrt{2}d}{\sigma_{n}} \right)}$ 00 $\frac{1}{4}{Q\left( \frac{3\sqrt{2}d}{\sigma_{n}} \right)}$ $\frac{1}{4}{Q\left( \frac{\sqrt{2}d}{\sigma_{n}^{2}} \right)}$ $\frac{1}{4}{Q\left( \frac{\sqrt{10}d}{\sigma_{n}} \right)}$ $\frac{1}{4}{Q\left( \frac{\sqrt{10}d}{\sigma_{n}} \right)}$ Average 1 ${\frac{1}{2}{Q\left( \frac{3\sqrt{2}d}{\sigma_{n}} \right)}} + {\frac{1}{2}{Q\left( \frac{\sqrt{2}d}{\sigma_{n}} \right)}}$ $Q\left( \frac{\sqrt{2}d}{\sigma_{n}} \right)$ ${\frac{1}{2}{Q\left( \frac{\sqrt{10}d}{\sigma_{n}} \right)}} + {\frac{1}{2}{Q\left( \frac{\sqrt{2}d}{\sigma_{n}} \right)}}$ ${\frac{1}{2}{Q\left( \frac{\sqrt{10}d}{\sigma_{n}} \right)}} + {\frac{1}{2}{Q\left( \frac{\sqrt{2}d}{\sigma_{n}} \right)}}$ Average 2 ${{\frac{1}{4}{Q\left( \frac{3\sqrt{2}d}{\sigma_{n}} \right)}} + {\frac{3}{4}{Q\left( \frac{\sqrt{2}d}{\sigma_{n}} \right)}}} \cong {\frac{3}{4}{Q\left( \frac{\sqrt{2}d}{\sigma_{n}} \right)}}$ ${{\frac{1}{2}{Q\left( \frac{\sqrt{10}d}{\sigma_{n}} \right)}} + {\frac{1}{2}{Q\left( \frac{\sqrt{2}d}{\sigma_{n}} \right)}}} \cong {\frac{1}{2}{Q\left( \frac{\sqrt{2}d}{\sigma_{n}} \right)}}$

The difference of bit error probability between MSB and LSB is decreased via swapping operation which is a significant improvement in terms of the average probability of bit error.

Through our proposed design criterion for horizontal bitwise mapping, we can find the optimal mapping sets for both 16-QAM and 64-QAM as shown in Table 5.

TABLE 5 # of trans. Mapping set # of trans. Mapping set 0 i₁q₁i₂q₂ 0 i₁q₁i₂q₂i₃q₃ 1 i₂q₂i₁q₁ 1 i₃q₃ī₂ q ₂i₁q₁ 2 i₁q₁ī₂ q ₂ 2 i₂q₂ī₃q₃i₁q₁ 3 i₂q₂ī₁ q ₁ 3 i₁q₁ī₂q₂i₃q₃ (a) 16 QAM (b) 64 QAM

The types of component bits of 64-QAM are divided into three sorts which are MSB, SB and LSB, respectively. The corresponding decision boundary on each component bit is introduced in FIG. 16. FIG. 16 shows decision boundary of 8-PAM. The difference of bit error probability between component bits within 64-QAM is larger than that of 16-QAM. Also, the difference of bit error probability between candidate component bits is significantly increased. For the first retransmission, averaging the difference of both the probability of bit error between MSB and LSB and the probability of bit error between candidate component bits corresponding to SB through swapping and inversion, respectively, is the best choice in terms of average bit error probability. For the second retransmission, employing second swapping operation between MSB and LSB after first swapping operation between MSB and SB with inversion operation of MSB has the best average bit error probability. Finally, in case of third retransmission, we find the optimal horizontal mapping that employs the inversion for SB through Equation 26.

Hereinafter, vertical bitwise mapping is disclosed.

HARQ schemes for multiple antenna system can exploit a spatial diversity. But, HARQ schemes such as Space Time Coding (STC)-HARQ are hard to obtain enough gain by spatial diversity due to time latency between retransmissions. Furthermore, because spatial diversity is only related to the increase of received SNR, the difference of reliability between component bits cannot be decreased by spatial selective diversity. It even aggravates the reliability difference.

In order to analyze the effect of vertical bitwise mapping, the channel model for 2 by 2 MIMO systems is defined by h_(a,b). The indexes of transmit antenna and retransmission are denoted by a and b, respectively. Under an assumption that the selectivity of time channel is small during retransmission, h_(a,b) can be represented by Equation 27.

MathFigure 27

h _(a,i) ≅h _(a,i+1) =h _(a)  [Math.27]

In addition to the case of time selective channel, low spatial selective channel can be also represented by Equation 28.

MathFigure 28

h _(a,i) ≅h _(a+1,i) =h _(a)  [Math.28]

Let us assume the low mobility channel under a rich scattering environment. In addition, 4-PAM for simpler analysis of vertical bitwise mapping is considered. The probability of error for a conventional HARQ-CC can be formulated as shown:

$\begin{matrix} {{MathFigure}\mspace{14mu} 29} & \; \\ {{\frac{1}{8}{Q\left( \frac{\sqrt{2{h_{1,i}}^{2}d}}{\sigma_{n}} \right)}} + {\frac{1}{4}{Q\left( \frac{\sqrt{2{h_{1,i}}^{2}d}}{\sigma_{n}} \right)}} + {\frac{1}{8}{Q\left( \frac{\sqrt{2{h_{2,i}}^{2}d}}{\sigma_{n}} \right)}} + {\frac{1}{4}{Q\left( \frac{\sqrt{2{h_{2,i}}^{2}d}}{\sigma_{n}} \right)}}} & \left\lbrack {{Math}.\mspace{14mu} 29} \right\rbrack \end{matrix}$

where i indicates the index of retransmission. We know that there is no additional gain extracted from either time or spatial channel diversity. Considering MSB shuffling between two transmit antennas during retransmission, the average probability of bit error is given as shown below.

$\begin{matrix} {{MathFigure}\mspace{14mu} 30} & \; \\ {{\frac{1}{4}{Q\left( \frac{\sqrt{{h_{1,i}}^{2} + {{h_{2,i}}^{2}d}}}{\sigma_{n}} \right)}} + {\frac{1}{2}{Q\left( \frac{\sqrt{2{h_{1,i}}^{2}d}}{\sigma_{n}} \right)}}} & \left\lbrack {{Math}.\mspace{14mu} 30} \right\rbrack \end{matrix}$

In case of LSB shuffling, the average probability of bit error is calculated as shown below.

$\begin{matrix} {{MathFigure}\mspace{14mu} 31} & \; \\ {{\frac{1}{4}{Q\left( \frac{\sqrt{2{h_{1,i}}^{2}d}}{\sigma_{n}} \right)}} + {\frac{1}{2}{Q\left( \frac{\sqrt{{h_{1,i}}^{2} + {{h_{2,i}}^{2}d}}}{\sigma_{n}} \right)}}} & \left\lbrack {{Math}.\mspace{14mu} 31} \right\rbrack \end{matrix}$

Finally, symbol shuffling is expressed as shown below.

$\begin{matrix} {{MathFigure}\mspace{14mu} 32} & \; \\ {{\frac{1}{4}{Q\left( \frac{\sqrt{{h_{1,i}}^{2} + {{h_{2,i}}^{2}d}}}{\sigma_{n}} \right)}} + {\frac{1}{2}{Q\left( \frac{\sqrt{{h_{1,i}}^{2} + {{h_{2,i}}^{2}d}}}{\sigma_{n}} \right)}}} & \left\lbrack {{Math}.\mspace{14mu} 32} \right\rbrack \end{matrix}$

In the above analysis, according to the number of bits that exchanged between antennas, the corresponding average bit error probability can be calculated by approach-based Q-function as same as for horizontal bitwise mapping. Through this analysis, we can know that MSB shuffling is the lowest average bit error probability. The following one of low bit error probability is LSB shuffling. The average bit error probability of symbol level shuffling is lower than that of conventional HARQ-CC, but is higher than that of 1-bit shuffling. In the similar way, we can compare with all cases of bit error probability derived by vertical mapping under a channel model and find a set of optimal vertical mapping.

Considering both horizontal bitwise mapping and vertical bitwise mapping, the proposed design criterion for the optimal bitwise mapping is represented as shown:

$\begin{matrix} {{MathFigure}\mspace{14mu} 33} & \; \\ {\left\{ {\gamma_{m},\mu_{m}} \right\} = {\arg \underset{\gamma_{i} \in \gamma}{\underset{\mu_{i} \in \mu}{Max}}{\begin{matrix} {{\frac{1}{2^{{L/2} - 1}L}{\sum\limits_{k}^{2^{L/2}}\; {\sum\limits_{j}^{L/2}\; {Q\left( \sqrt{\frac{1}{2\sigma_{n}^{2}}{\sum\limits_{i = 1}^{m}\; {{{\gamma_{i}\left( h_{j,k} \right)}}^{2}{{\mu_{i}\left( d_{j,k} \right)}}^{2}}}} \right)}}}} -} \\ {\frac{1}{2^{{L/2} - 1}L}{\sum\limits_{k}^{2^{L/2}}\; {\sum\limits_{j}^{L/2}\; {Q\left( \sqrt{\frac{1}{2\sigma_{n}^{2}}{\sum\limits_{i = 1}^{m - 1}\; {{{\gamma_{i}\left( h_{j,k} \right)}}^{2}{{\mu_{i}\left( d_{j,k} \right)}}^{2}}}} \right)}}}} \end{matrix}}}} & \left\lbrack {{Math}.\mspace{14mu} 33} \right\rbrack \end{matrix}$

where γ_(i) denotes a vertical mapping at i-th retransmission. γ_(i) subject to a set of vertical mapping γ. As well as the way to reduce the size of set in Equation 26, we can reduce the size of set of candidate vertical mapping in course of finding the set of optimal vertical mapping as excluding the cases with low priority. From our design criterion for a combination of horizontal bitwise mapping and vertical bitwise mapping, Table 6 shows the set of optimal mapping.

TABLE 6 # of re- # of re- transmission Optimal mapping transmission Optimal mapping 0 i_(1,1)q_(1,1)i_(1,2)q_(1,2) 0 i_(1,1)q_(1,1)i_(1,2)q_(1,2)i_(1,3)q_(1,3) i_(2,1)q_(2,1)i_(2,2)q_(2,2) i_(2,1)q_(2,1)i_(2,2)q_(2,2)i_(2,3)q_(2,3) 1 i_(1,2)q_(1,2)i_(2,1)q_(2,1) 1 i_(2,3)q_(1,3)ī_(2,2) q _(1,2)i_(2,1)q_(1,1) i_(2,2)q_(2,2)i_(1,1)q_(1,1) i_(1,3)q_(2,3)ī_(1,2) q _(2,2)i_(1,1)q_(2,1) 2 i_(1,1)q_(1,1)ī_(2,2) q _(2,2) 2 i_(2,2)q_(2,2)ī_(2,3) q _(2,3)i_(2,1)q_(1,1) i_(2,1)q_(2,1)ī_(1,2) q _(1,2) i_(1,2)q_(1,2)ī_(1,3) q _(1,3)i_(1,1)q_(2,1) 3 i_(1,2)q_(1,2)ī_(2,1) q _(2,1) 3 i_(2,1)q_(1,1)ī_(2,2) q _(1,2)i_(2,3)q_(1,3) i_(2,2)q_(2,2)ī_(1,1) q _(1,1) i_(1,1)q_(2,1)ī_(1,2) q _(2,2)i_(1,3)q_(2,3) (a) 16-QAM (b) 64-QAM

The proposed bitwise mapping is simulated for 16-QAM with Convolution Turbo Coding (CTC) of rate ½ under a 2 by 2 MIMO channel. The most parameter for performance evaluation follows the standard specification of IEEE 802.16e. For channel model, the typical urban (TU) is considered as frequency selective MIMO channel. Also, time latency between retransmissions is 5 ms and the number of retransmissions is considered up to 3. Full usage of the sub-channels (FUSC) which is one of distributed resource allocation is considered to obtain the coding gain on frequency domain as sufficient as possible. The type of receiver is based on linear Minimum Mean Square Error (MMSE) equalizer. For 10 MHz bandwidth, the number of subcarriers is 1024.

FIG. 17 shows Frame Error Rate (FER) performances according to the number of bits shuffled between transmit antennas. Simulation is performed on all possible combinations which can be exchanged between transmit antennas by component bits within 64-QAM. Based on the analysis for the vertical bit mapping, we can know that probability of average bit error for SB shuffling between transmit antennas has the best performance at 1% FER. The result of simulation also corresponds to the analysis of probability of bit error-based Q-function.

FIG. 18 compares the FER performance of optimal bitwise mapping and conventional mapping (i.e. HARQ-CC). As opposed to the typical scheme in which the combining gain is decreased, as the number of retransmissions is increased, the proposed optimal bitwise mapping always obtains the combining gain of almost about 2 dB regardless of the number of retransmissions. In case of third retransmission, the proposed scheme provides about 8 dB gain over the conventional mapping. This performance is due to the coding gain which can be obtained by diversity gain extracted from time selectivity and spatial selectivity through proposed bitwise mapping.

III. Constellation Rearrangement for IR Resource

In case of circular buffer-based HARQ-incremental redundancy (IR), a redundancy version is transmitted in every retransmission. Coded bits, which are the same as those previously transmitted, can be transmitted when the redundancy version is transmitted. In this case, constellation rearrangement can be performed for these coded bits. A swapping operation and/or an inversion operation can be performed for constellation rearrangement.

FIG. 19 shows HARQ-IR using a swapping operation. A cyclic buffer 400 stores a code block. An initial data block 410 is a part of the code block and is initially transmitted. According to a retransmission request of the initial data block 410, in 1st retransmission, a 1st retransmission block 420 subsequent to the initial data block 410 in the cyclic buffer 400 is transmitted. In 2nd retransmission, a 2nd retransmission block 430 consists of a non-wraparound block 430 a and a wraparound block 430 b according to a characteristic of the cyclic buffer 400. The wraparound block 430 b denotes a block which has previously been transmitted (or retransmitted) one time when HARQ is performed. The wraparound block 430 b included in the 2nd retransmission block 430 is wrapped around with a part of the initial data block 410. By performing the swapping operation for constellation rearrangement on the wraparound block 430 b, reliability differences for all bits are averaged. In case of 16-QAM, bit reliabilities of MSB and LSB can be swapped and then averaged. In case of 64?QAM, bit reliabilities of MSB, SB, and LSB can be swapped and then averaged. In 3rd retransmission, since all parts of a 3rd retransmission block 440 belong to a data block which has previously been transmitted one time, the 3rd retransmission block 440 is a wraparound block and thus constellation rearrangement is performed thereon. In 4th retransmission, since all parts of a 4th retransmission block 450 belong to a data block which has previously been transmitted one time, the 4th retransmission block 450 is a wraparound block and thus constellation rearrangement is performed thereon.

The size of cyclic buffer and the size of redundancy version transmitted according to the maximum number of retransmissions are for exemplary purposes only, and thus the technical features of the present invention are not limited thereto. In addition, the swapping operation for constellation rearrangement is also for exemplary purposes only, and thus an inversion operation, a spatial operation, and a combination of the two operations can also be performed.

FIG. 20 shows HARQ-IR using a swapping operation and an inversion operation. A cyclic buffer 500 stores a code block. An initial data block 510 is a part of the code block and is initially transmitted. According to a retransmission request of the initial data block 510, in 1st retransmission, a 1st retransmission block 520 subsequent to the initial data block 510 in the cyclic buffer 500 is transmitted. The 1st retransmission block 520 consists of a non-wraparound block 520 a and a wraparound block 520 b according to a characteristic of the cyclic buffer 500. This is because the wraparound block 520 b included in the 1st retransmission block 520 is wrapped around with a part of the initial data block 510. The swapping operation for constellation rearrangement is performed on the wraparound block 520 b. In 2nd retransmission, a 2nd retransmission block 530 consists of one-time wraparound block 530 a and a two-time wraparound block 530 b. The one-time wrapped block 530 a denotes a block which has previously been transmitted one time in a wraparound manner. The two-time wrapped block 530 b denotes a block which has previously been transmitted two times in a wraparound manner. A constellation rearrangement pattern may vary according to the number of wraparounds. Herein, the one-time wrapped block 530 a uses the swapping operation, and the two-time wrapped block 530 b uses the inversion operation. All parts of a 3rd retransmission block 540 are 2-time wraparound blocks, and the inversion operation is performed thereon. The constellation rearrangement pattern depends on the number of wraparounds and thus may vary according to the number of wraparounds. However, the technical features of the present invention are not limited to the number of wraparounds or the constellation rearrangement pattern.

Tables 7 and 8 below show bitwise mapping for each modulation scheme according to the number of wraparounds in a process of retransmitting coded bits in a cyclic buffer. 16-QAM is used in Table 7. 64-QAM is used in Table 8.

TABLE 7 # of wraparound Mapping rule 0 i₁ q₁ i₂ q₂ 1 i₂ q₂ i₁ q₁ 2 i₁ q₁ ī₂ q ₂ 3 i₂ q₂ ī₁ q ₁

TABLE 8 # of wraparound Mapping rule 0 i₁ q₁ i₂ q₂ i₃ q₃ 1 i₃ q₃ ī₂ q ₂ i₁ q₁ 2 i₂ q₂ ī₃ q ₃ i₁ q₁ 3 i₁ q₁ ī₂ q ₂ i₃ q₃

A single antenna system is considered in the HARQ-IR. For the HARQ-IR of the MIMO system, spatial bitwise rearrangement can also be performed to obtain spatial diversity. Constellation rearrangement is performed between transmit antennas whenever coded bits in a cyclic buffer are transmitted in a wraparound manner in a retransmission process. For example, when retransmission is first made in a wraparound manner, a swapping operation is performed on the signal constellation and also spatial bitwise mapping is performed between antennas. When retransmission is secondly made in a wraparound manner, an inversion operation is performed on the signal constellation and also spatial bitwise mapping is performed between antennas.

Tables 9 and 10 below show bitwise mapping for each modulation scheme according to the number of wraparounds in a process of retransmitting coded bits in a cyclic buffer in a multiple antenna system having two transmit antennas. 16-QAM is used in Table 9. 64-QAM is used in Table 10. In i_(a,b) and q_(a,b), ‘a’ denotes an antenna index, and ‘b’ denotes an index indicating coordinates on the constellation.

TABLE 9 # of wraparound Mapping rule 0 i_(1,1)q_(1,1)i_(1,2)q_(1,2) i_(2,1)q_(2,1)i_(2,2)q_(2,2) 1 i_(1,2)q_(1,2)i_(2,1)q_(2,1) i_(2,2)q_(2,2)i_(1,1)q_(1,1) 2 i_(1,1)q_(1,1)ī_(2,2) q _(2,2) i_(2,1)q_(2,1)ī_(1,2) q _(1,2) 3 i_(1,2)q_(1,2)ī_(2,1) q _(2,1) i_(2,2)q_(2,2)ī_(1,1) q _(1,1)

TABLE 10 # of wraparound Mapping rule 0 i_(1,1)q_(1,1)i_(1,2)q_(1,2)i_(1,3)q_(1,3) i_(2,1)q_(2,1)i_(2,2)q_(2,2)i_(2,3)q_(2,3) 1 i_(2,3)q_(1,2)ī_(2,2) q _(1,2)i_(2,1)q_(1,1) i_(1,3)q_(2,3)ī_(1,2) q _(2,2)i_(1,1)q_(2,1) 2 i_(2,2)q_(2,2)ī_(2,3) q _(2,3)i_(2,1)q_(1,1) i_(1,2)q_(1,2)ī_(1,3) q _(1,3)i_(1,1)q_(2,1) 3 i_(2,1)q_(1,1)ī_(2,2) q _(1,2)i_(2,3)q_(1,3) i_(1,1)q_(2,1)ī_(1,2) q _(2,2)i_(1,3)q_(2,3)

To obtain frequency diversity, when symbols are mapped to subcarriers in a retransmission process, mapping can be extended towards noncontiguous subcarriers to perform bitwise rearrangement on a signal constellation. Frequency diversity can be additionally obtained by extending mapping towards the noncontiguous subcarriers. In the mapping from the symbols to the subcarriers, an additional HARQ gain can be ensured by applying constellation rearrangement to obtain an average reliability of coded bits.

FIG. 21 is a block diagram showing an apparatus for wireless communication according to an embodiment of the present invention. An apparatus 50 for wireless communication may be a part of a UE. The apparatus 50 for wireless communication includes a processor 51, a memory 52, a radio frequency (RF) unit 53, a display unit 54, and a user interface unit 55. The RF unit 53 is coupled to the processor 51 and transmits and/or receives radio signals. The memory 52 is coupled to the processor 51 and stores an operating system, applications, and general files. The display unit 54 displays a variety of information of the UE and may use a well-known element such as a liquid crystal display (LCD), an organic light emitting diode (OLED), etc. The user interface unit 55 can be configured with a combination of well-known user interfaces such as a keypad, a touch screen, etc. With an implemented physical layer, the processor 51 supports the aforementioned constellation rearrangement and HARQ. The proposed method can be implemented by the processor 51.

All functions described above may be performed by a processor such as a microprocessor, a controller, a microcontroller, and an application specific integrated circuit (ASIC) according to software or program code for performing the functions. The program code may be designed, developed, and implemented on the basis of the descriptions of the present invention, and this is well known to those skilled in the art.

While the present invention has been particularly shown and described with reference to exemplary embodiments thereof, it will be understood by those skilled in the art that various changes in form and details may be made therein without departing from the spirit and scope of the invention as defined by the appended claims. The exemplary embodiments should be considered in descriptive sense only and not for purposes of limitation. Therefore, the scope of the invention is defined not by the detailed description of the invention but by the appended claims, and all differences within the scope will be construed as being included in the present invention. 

1. A data retransmission method using hybrid automatic repeat request (HARQ), comprising: transmitting a data block; receiving a retransmission request signal for the data block; generating a retransmission block by performing swapping or inversion between bits constituting the data block according to the retransmission request signal; and transmitting the retransmission block.
 2. The method of claim 1, wherein the swapping between bits is performed by swapping a bit having a high bit reliability and a bit having a low bit reliability in previous transmission.
 3. The method of claim 1, wherein the data block is a transmission block having a matrix format and consisting of rows whose number is equal to the number of transmit antennas, each row of the transmission block is transmitted through one transmit antenna, bits constituting one row represent bits of one data symbol, and the data symbol is modulated using at least two modulation schemes.
 4. The method of claim 3, wherein the swapping between bits is performed by cyclically shifting bits constituting each column of the transmission block, and each column has a different cyclic shift amount.
 5. The method of claim 1, wherein the retransmission block is mapped according to the number of retransmissions when using a 16-ary quadrature amplitude modulation (QAM) or 64-QAM scheme as expressed in the following table: # of trans. Mapping set # of trans. Mapping set 0 i₁q₁i₂q₂ 0 i₁q₁i₂q₂i₃q₃ 1 i₂q₂i₁q₁ 1 i₃q₃ī₂ q ₂i₁q₁ 2 i₁q₁ī₂ q ₂ 2 i₂q₂ī₃q₃i₁q₁ 3 i₂q₂ī₁ q ₁ 3 i₁q₁ī₂q₂i₃q₃ (a) 16 QAM (b) 64 QAM

where, in i_(a) and q_(a) of the table, ‘a’ denotes an index indicating coordinates on a constellation.
 6. The method of claim 1, wherein the retransmission block is mapped according to the number of retransmissions when using a 16-QAM or 64-QAM scheme as expressed in the following table: # of re- # of re- transmission Optimal mapping transmission Optimal mapping 0 i_(1,1)q_(1,1)i_(1,2)q_(1,2) 0 i_(1,1)q_(1,1)i_(1,2)q_(1,2)i_(1,3)q_(1,3) i_(2,1)q_(2,1)i_(2,2)q_(2,2) i_(2,1)q_(2,1)i_(2,2)q_(2,2)i_(2,3)q_(2,3) 1 i_(1,2)q_(1,2)i_(2,1)q_(2,1) 1 i_(2,3)q_(1,3)ī_(2,2) q _(1,2)i_(2,1)q_(1,1) i_(2,2)q_(2,2)i_(1,1)q_(1,1) i_(1,3)q_(2,3)ī_(1,2) q _(2,2)i_(1,1)q_(2,1) 2 i_(1,1)q_(1,1)ī_(2,2) q _(2,2) 2 i_(2,2)q_(2,2)ī_(2,3) q _(2,3)i_(2,1)q_(1,1) i_(2,1)q_(2,1)ī_(1,2) q _(1,2) i_(1,2)q_(1,2)ī_(1,3) q _(1,3)i_(1,1)q_(2,1) 3 i_(1,2)q_(1,2)ī_(2,1) q _(2,1) 3 i_(2,1)q_(1,1)ī_(2,2) q _(1,2)i_(2,3)q_(1,3) i_(2,2)q_(2,2)ī_(1,1) q _(1,1) i_(1,1)q_(2,1)ī_(1,2) q _(2,2)i_(1,3)q_(2,3) (a) 16-QAM (b) 64-QAM

where, in i_(a,b) and q_(a,b) of the table, ‘a’ denotes an antenna index and ‘b’ denotes an index indicating coordinates on a constellation.
 7. A data retransmission method in a multiple-antenna system, comprising: transmitting a first data symbol modulated with a first modulation scheme through a first transmit antenna and transmitting a second data symbol modulated with a second modulation scheme through a second transmit antenna; generating a new first data symbol and a new second data symbol by performing swapping or inversion on bits constituting the first data symbol and bits constituting the second data symbol; and transmitting the new first data symbol through the first transmit antenna and transmitting the new second data symbol through the second transmit antenna.
 8. The method of claim 7, further comprising receiving a retransmission request signal for the first data symbol and the second data symbol.
 9. An apparatus for wireless communication, comprising: a radio frequency (RF) unit for transmitting and receiving a radio signal; and a processor coupled with the RF unit, for performing HARQ-increment redundancy (IR), and configured to: sequentially configure data blocks in a cyclic buffer of the HARQ-IR according to a retransmission request, and transmit the data blocks, wherein if all or some parts of the data block are wrapped around with a previously transmitted data block, constellation rearrangement is performed on the wraparound parts.
 10. The apparatus of claim 9, wherein the processor is configured to use 16-QAM as a modulation scheme, and perform constellation rearrangement according to the number of wraparounds as described in the following table: # of wraparound Mapping rule 0 i₁ q₁ i₂ q₂ 1 i₂ q₂ i₁ q₁ 2 i₁ q₁ ī₂ q ₂ 3 i₂ q₂ ī₁ q ₁

where, in i_(a) and q_(a) of the table, ‘a’ denotes an index indicating coordinates on a constellation.
 11. The apparatus of claim 9, wherein the processor is configured to use 16-QAM as a modulation scheme, and perform constellation rearrangement according to the number of wraparounds as described in the following table: # of wraparound Mapping rule 0 i_(1,1)q_(1,1)i_(1,2)q_(1,2) i_(2,1)q_(2,1)i_(2,2)q_(2,2) 1 i_(1,2)q_(1,2)i_(2,1)q_(2,1) i_(2,2)q_(2,2)i_(1,1)q_(1,1) 2 i_(1,1)q_(1,1)ī_(2,2) q _(2,2) i_(2,1)q_(2,1)ī_(1,2) q _(1,2) 3 i_(1,2)q_(1,2)ī_(2,1) q _(2,1) i_(2,2)q_(2,2)ī_(1,1) q _(1,1)

where, in i_(a,b) and q_(a,b) of the table, ‘a’ denotes an antenna index and ‘b’ denotes an index indicating coordinates on a constellation.
 12. The apparatus of claim 9, wherein the processor is configured to use 64-QAM as a modulation scheme, and perform constellation rearrangement according to the number of wraparounds as described in the following table: # of wraparound Mapping rule 0 i₁q₁i₂q₂i₃q₃ 1 i₃q₃ī₂ q ₂i₁q₁ 2 i₂q₂ī₃ q ₃i₁q₁ 3 i₁q₁ī₂ q ₂i₃q₃

where, in ia and qa of the table, ‘a’ denotes an index indicating coordinates on a constellation.
 13. The apparatus of claim 9, wherein the processor is configured to use 64-QAM as a modulation scheme, and perform constellation rearrangement according to the number of wraparounds as described in the following table: # of wraparound Mapping rule 0 i_(1,1)q_(1,1)i_(1,2)q_(1,2)i_(1,3)q_(1,3) i_(2,1)q_(2,1)i_(2,2)q_(2,2)i_(2,3)q_(2,3) 1 i_(2,3)q_(1,3)ī_(2,2) q _(1,2)i_(2,1)q_(1,1) i_(1,3)q_(2,3)ī_(1,2) q _(2,2)i_(1,1)q_(2,1) 2 i_(2,2)q_(2,2)ī_(2,3) q _(2,3)i_(2,1)q_(1,1) i_(1,2)q_(1,2)ī_(1,3) q _(1,3)i_(1,1)q_(2,1) 3 i_(2,1)q_(1,1)ī_(2,2) q _(1,2)i_(2,3)q_(1,3) i_(1,1)q_(2,1)ī_(1,2) q _(2,2)i_(1,3)q_(2,3)

where, in i_(a,b) and q_(a,b) of the table, ‘a’ denotes an antenna index and ‘b’ denotes an index indicating coordinates on a constellation. 